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Abstract: We use WMAP5 and other cosmological data to constrain model parameters 
in quintessence cosmologies, focusing also on their shift when we allow for non-vanishing 
neutrino masses. The Ratra-Peebles (RP) and SUGRA potentials are used here, as ex- 
amples of slowly or fastly varying state parameter w{a). Both potentials depend on an 
energy scale A. Here we confirm the results of previous analysis with WMAP3 data on the 
upper limits on A, which turn out to be rather small (down to ~ 10~^ in RP cosmologies 
and ~ 10~^ for SUGRA). Our constraints on A are not heavily affected by the inclusion 
of neutrino mass as a free parameter. On the contrary, when the neutrino mass degree of 
freedom is opened, significant shifts in the best-fit values of other parameters occur. 
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1. Introduction 

Recent data require an accelerated expansion of the Universe. High-redshift supernovae Q, 
CMB [§] and deep samples essentially agree with requiring a density parameter for non- 
relativistic matter ~ 0.27, while the total density parameter f^o ~ 1- A fluid, dubbed 
Dark Energy (DE), filling the gap in the cosmological energy balance, can account for 
cosmic acceleration if its state parameter w = Pde/Pde is close to -1. 

In the minimal model, dubbed ACDM (or cosmic concordance) cosmology, DE has a 
state parameter 11! = —1 as though being a false vacuum, and this is equivalent to intro- 
ducing Einstein's cosmological constant. All data available up to now can be accomodated 
in this cosmology. Still there are several problems related to the nature of the DE. Infact, 
two critical issues plague such model: the fine tuning of the present DE energy density to 
the vacuum expectation value of quantum field theories, and the coincidence between DE 
domination and the advent of the non-linear growth of matter fluctuations. 

In the attempt to approach a better understanding of DE nature, a varying state 
parameter w{a) has been tentatively allowed. Here two patterns can be followed: either 
phaenomenological fitting expressions for w{a) can be introduced or, if we assume that DE 
is a self-interacting scalar field (/> tracking potential expressions V{(l)) (5|, depending on 
suitable parameters, can be tested. This latter approach has the advantage of allowing for 
larger energy scales of DE, thus easing the fine tuning problem. 

Analysis of such dynamical Dark Energy (dDE) scenario were previously performed 
when WMAP3 data became available Q. The first aim of this work is to inspect how far 
their output are modified when WMAP5 data are considered, in analogy to the modifica- 
tions occurring for ACDM cosmologies. As in j^, we shall do so by using RP and SUGRA 
potentials (see below). 

Here we shall however further focus on the impact on parameters arising from allowing 
for a neutrino mass. As a matter of fact, cosmology is sensitive to u masses (see Q, [|8|) 
even at a level below current limits set by /3-decay experiments |^], so much that the very 
opening of the degree of freedom of neutrino mass, causes appreciable shifts in the best-fit 
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values of some cosmological parameters. It should be outlined that non-vanishing neutrino 
masses are required by flavor mixing results [jl^ . Even though such results set quite a low 
limit to the sum My of neutrino masses, it would be therefore unappropriate to neglect the 
neutrino mass option. 

This fact was recently exploited by several researchers to show that a larger coupling 
between CDM and DE is compatible with data [0, 12, ^] when allowing for a non-zero 



neutrino mass. Here, however, no such coupling shall be considered. 

Let us finally outline that further alternatives to a cosmological constant or a scalar 
field have been considered. Among them, modifications to Hilbert-Einstein Lagrangian 
density, adding f{R) terms {R: Ricci scalar), or the hypothesis that DE arises as a back 



reaction to the development of inhomogeneities |14|. Such options make the need of an 
extra substance superfluous, but give rise to other theoretical and observational problems. 
We shall not deal with them in this work. 

The plane of the paper is as follow. In Sec. ^ we describe the models. Data and 
methods used are reviewed in Sec. |3|, while in Sec. ^ we present results and discuss our 
flndings. 

2. Models 

In a dDE scenario and in the reference frame where the metric is ds^ = a^(T)[dT^ — di'^] 
(r : conformal time; di : spatial line element), the energy density and pressure of the held 
(f) read 

Pde = Pk + V{(l)) , Pde = Pk- V{<P) , with pk = {(t>'f/2a^ , (2.1) 

prime indicating differentiation in respect to r. Accordingly, \i pk ^ V , the DE state 
parameter approaches +1 {stiff matter) so that DE energy density rapidy dilutes during 
expansion [p oa a~^). In the opposite case V ^ pk, the state parameter approaches — 1 
and DE is suitable to explain the observed cosmic acceleration. 
The equations 

p',, + 3-(P,e + Pde)=0, (2.2) 



a 

and ^ 

(P" + 2%' + a^V,^=0 (2.3) 

are then clearly equivalent. By integrating ( |2.3[) together with the Friedmann equations, 
we obtain the time evolution of (j) and, thence, of pk, V{(j)), pde, Pde and, flnally, t(;(a). 

Quite in general, the solution of a differential equation depends on initial conditions. 
Tracker potentials, however, have attractor solutions (tracking solutions) on which they 
converge (almost) independently of the initial conditions. 

Among tracker potentials, we consider here: 

V{<p) = A"+V</'" RP |ll (2.4) 

y(0) = (A°+V(/)")exp(47r</.Vm^) SUGRA ||]. (2.5) 
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Figure 1: Comparison between the variation of the state parameter for RP and SUGRA potentials 
and the usual parametrization w(a) — Wo + {l+a)wa- Here, A is 100 GeV for SUGRA and 10~* GeV 
for RP, while {wo, Wa) — (-0.8908, 0.525). These values are selected to yield a behavior close to 
SUGRA, at low z. Even renouncing to a full coincidence at z = 0, the fast variablility of w{a) in 
SUGRA cannot be met by any polinomial w{a). Notice also that, for a RP potential, w{a) varies 
more slowly. The plot is for h = 0.7, n,, = 0.046, = 0.209. 

The RP potential yields a slowly varying w{a) state parameter, while the variation of w{a) 
for SUGRA potential is faster, as shown in Figure [l|. In both cases, for any choice of A 
and a these potentials yield a precise DE density parameter i^^e- Here we use A and Qde 
as free parameters in flat cosmologies; the related a value is then suitably fixed. 

Let us then remind the constraints on RP and SUGRA potentials obtained when using 
WMAP3 data |^]. They confirmed that a RP cosmology requires an energy scale A below 
~ 10~'^GeV (and then a ~ 0.5); this is coherent with the fact that w^a) is then steadily 
below -0.85 (~ 95% C.L.) as required from observations for a constant w On the 

contrary, in the SUGRA case, while their best fit still leads to low A's, they found that 
a value A ~ 100 GeV is within ~ 2 it's from the best-fitting model; such scale is then 
compatible with the ones of SUSY or EW transition. 

One significant result of our analysis is that these upper limits are badly lowered when 
using WMAP5 data. 



3. Methods 



For this work we use a modified version of the cosmological Boltzmann code CAMB |18, 19 
to calculate theoretical CMB and matter power spectra. In our modified code, a precursor 
program is called before running CAMB, evolving the background solutions in order to 
determine iteratively the value of a consistent with the assigned ilde and A. The scalar 
field (/) evolves according to the Klein-Gordon equation (p.3|). Its initial value (j)o is obtained 
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considering that pde,o ~ ^{4'o)- In the code we also took into account spatial fluctuations 
of the DE field. 

Our code is embedded into the publicly available Monte Carlo Markov Chain (MCMC) 
engine CosmoMC |2^]. The following set of free parameters are used: {ub, ^c, d, To, ng, 
ln{W^^As), A, M^}. Here LOb,c = ^b,ch'^ are the reduced density parameters for baryonic and 
cold dark matter, respectively, h being the present value of the Hubble constant in units of 
100 km/s/Mpc; 9 is then the ratio of the sound horizon to the angular diameter distance at 
recombination, Tq is the optical depth, Ug is the scalar spectral index, Ag is the primordial 
amplitude of scalar fluctuations at a scale k = 0.05Mpc~^ and A = logio(A/GeV). In 
addition we marginalize over the SZ amplitude. We assume three massive neutrino species 
with degenerate masses. 

In the MCMC runs we combined the five year data from the WMAP measurements of 
the CMB radiation (WMAP5) with the galaxy power spectrum from the 2dF survey pT|| , 
supernovae type la data from the SNLS survey [p^], and added the priors: h = 0.72 it 0.08 



(from the HST key project H), = 0.022 ± 0.002 (from BBN results ||, |2|, 
4. Results and discussion 

In Table |l] we report best fit parameter values and 1-a errors for different models. In the 
case of neutrino masses, whose sum Mi, is not determined by data, we give its (95 % C.L.) 
upper limit. Also for the scales A = log(A/GeV), in the RP and SUGRA potentials, an 
upper limit is given. 

Let us first comment on the ACDM and w = const, column, which was added here for 
the sake of comparison. In both of them neutrino masses are allowed. As expected, the 
values we find for all parameters are strictly coeherent with those found by the WMAP 
team when combining the WMAP5 data with measurements of the position of the baryonic 
acoustic oscillation (BAO) peak in the galaxy power spectrum and SNLS data. The main 
discrepancy concerns error bars, being systematically slightly greater here than in the 
outputs from the WMAP team. 

The reason for this can be easily understood. Deep sample data as 2dF and SDSS are 
substantially coherent, as far as the BAO location is concerned. Some discrepancy between 
them is however present in the spectral details. Also because of that, the WMAP5 fit takes 
into account the BAO position only. On the contrary, here we took into account the whole 
information contained in the 2dF spectral data. In addition we have used slightly different 
priors than the WMAP team did, so that the (small) discrepancies found are fully justified. 
In particular, our results for a constant w = —1 are largely consistent with the analysis of 
the WMAP team. 

Infact, for w = -Iwe find an upper limit of < 0.66eV, while WMAP5+SNLS+BA0 
yields Mi, < 0.61eV. However, when w ^ —1 models are considered, the upper limit ex- 
hibits a sharp increase, up to 0.94 eV. WMAP5+SNLS+BA0 outputs show only a ~ 10 % 
shift in the upper limit when passing from ACDM to w ^ —1. Here the shift is ~ 4 times 
greater. This is related to the widening of the w interval, as greater Mi^ values become 
compatible with data if we delve into the phantom regime |27|. Accordingly, as RP and 
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Parameter 



ACDM w = const. RP SUGRA 

/ M^^O M^ = / = JVU ^ 



2.258 2.247 2.278 2.271 2.278 2.274 

± 0.061 ± 0.062 ± 0.060 ± 0.060 ± 0.060 ± 0.060 

0.1098 0.1132 0.1051 0.1062 0.1043 0.1055 

± 0.0040 ± 0.0069 ± 0.0050 ± 0.0047 ± 0.0051 ± 0.0050 

, ^ 70.1 69.7 70.9 69.8 70.6 69.6 

(km/s/Mpc) ^2.2 ±1.7 ±2.0 ±1.9 ±2.0 



10^ LOb 



o 



0.087 0.085 0.089 0.090 0.089 0.091 

± 0.017 ± 0.017 ± 0.017 ± 0.018 ± 0.017 ± 0.018 



w -1 



-1.06 
± 0.10 



(eV) 
(95% C.L.) 
A 

(95% C.L.) 



n.. 



ln(10i%) 



< 0.66 < 0.94 - < 0.58 - < 0.59 

- - < -8.4 < -8.7 < -3.4 < -4.7 

0.962 0.958 0.968 0.967 0.970 0.967 

± 0.014 ± 0.015 ± 0.014 ± 0.014 ± 0.014 ± 0.014 

3.045 3.049 3.047 3.041 3.045 3.041 

± 0.040 ± 0.040 ± 0.041 ± 0.042 ± 0.041 ± 0.041 

0.270 0.280 0.254 0.266 0.255 0.266 

± 0.022 ± 0.027 ± 0.017 ± 0.021 ± 0.018 ± 0.020 

0.713 0.711 0.749 0.697 0.737 0.689 

± 0.056 ± 0.059 ± 0.039 ± 0.052 ± 0.046 ± 0.054 

-2 ln(£) 1407.25 1407.38 1407.35 1407.38 1407.44 1407.31 

Table 1: Constraints on ACDM and w =const models, with massive i^'s, and RP and SUGRA 
models, with and without massive z/'s. Limits obtained using the combination of data sets described 
in the text. The last row gives the likelihood of the best-fit sample for each model. 
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SUGRA models yield an effective equation of state w{a) > — 1, this leads to more stringent 
limits on My. 
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As far as RP and SUGRA models are concerned, the decrease of errors in respect to 
WMAP3 outputs in |^] has direct consequences on the hmits on the energy scale A. In the 
case Mjy = 0, treated in Q, an upper limit A < —7.7 (2.1) was found for RP (SUGRA) at 
95% C.L.. Such limits are now lowered to -8.4 for RP and to a severe -3.4 in the case of 
SUGRA. Limits are even more restrictive if M^, 7^ is allowed. 

The reason for this strong improvement in the limits on A can partly be understood 
by inspecting Figures § and|^. Here we show 2D contours for A versus lOc, erg, Hq and M,^, 
which are parameters that are strongly correlated with A. In the plots we compare the 
results from the models with M,^ = 0, to those with a free neutrino mass. There is, for 
example, strong anticorrelations between ag and A and lOc and A. As our analysis gives a 
preference of significantly larger ul>c and erg than the WMAP3 analysis in j^, we also find 
much tighter limits on A. Parts of the improved limits can also be attributed to our general 
decrease in error bars on the parameters. 




Figure 2: 2D contour plots with A vs Wc, erg, Hq and Mi, for models with a RP potential. Solid 
black lines are for models free M^, and dotted red lines show the contours for models with M^, = 



As far as the other parameters are concerned, all shifts in respect to ACDM are within 
1-a. It may be however worth considering the 2-cr upper limit for Ug. When using WMAP3 
data there was a discrepancy between the ACDM and the RP or SUGRA cases. In the 
former case Ug = 1 was excluded at the 95 % C.L.; on the contrary, = 1 kept consistent 
with data, at the same C.L., for RP and SUGRA. Using WMAP5 data the discrepancy 
disappears and = 1 is above the 95% C.L. for all models considered. This means that 
the related constraint on inflationary models is apparently much more independent from 
DE nature. 
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Figure 3: Same as Figure ^, but with a SUGRA potential. 

Another significant shift worth outlining concerns dg values. Let us remind that 

a| = A Jdk k'^+''^T\k)W\kRs) . (4.1) 

Here: T(k) is the transfer function, W{kR) is the Fourier transform of a (top-hat) filter, 
i?8 = S/i-^Mpc . 

Although in the absence of major shifts on Ag and n^, it is sufficient to open the degree 
of freedom of neutrino mass, to allow for significant shifts in T which, in turn, reflect onto 
a smaller best-fit cxg value. 

Also likelihood values for the best-fitting models do not exhibit major changes and 
cannot be used to discriminate among the cosmologies considered here. This is coherent 
with the general conclusion we draw, that current data are still unable favor any kind of 
DE nature. 

To summarize, we have studied dDE models with RP and SUGRA potentials, and 
used a compilation of cosmological data sets to constrain parameters of these models. We 
have also studied the effect of opening an extra degree of freedom, by allowing for non-zero 
neutrino masses. 

Our basic result is that the energy scale of the DE model can be constrained to A < —8.4 
(-3.4) if Mjy = and to A < —8.7 (-4.7) in presence of massive v, when using a RP (SUGRA) 
potential. These limits are significantly stronger than what was found in Q, where they 
used WMAP3 data and vanishing neutrino masses only. Some of the improvement can 
be attributed to the inclusion of neutrino masses, as non-zero neutrino masses tend to 
shrink the allowed parameter space for A (see Figures |2| and ^). However, the most of 
the improvement is caused by our use of more precise data, and that these data have a 
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preference of slightly different values of parameters like ag and lOc, which in turn drives A 
to smaller values. 
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